IISc-CTS-3/94 



hep-ph/94052151 



FROM KONDO MODEL AND STRONG COUPLING LATTICE QCD 
TO THE ISGUR-WISE FUNCTION^ 



(N 



APOORVA PATEL* 
CTS and SERC, Indian Institute of Science 
Bangalore-560012, India 



Q\ ■ ABSTRACT 



Isgur-Wise functions parametrise the leading behaviour of weak decay form factors 

of mesons and baryons containing a single heavy quark. The form factors for the 

quark mass operator are calculated in strong coupling lattice QCD, and Isgur-Wise 

functions extracted from them. Based on renormalisation group invariance of the 

operators involved, it is argued that the Isgur-Wise functions would be the same in 

the weak coupling continuum theory. 
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Quantum Chromodynamics with heavy quarks possesses spin-flavour symme- 
tries that become exact as the quark masses go to infinity. These symmetries 
give rise to relations amongst various matrix elements and form factors of hadrons 
. containing heavy quarks 1 . Such relations based on symmetry properties alone are 
genuine predictions of QCD, and do not suffer from the uncertainties of phenomeno- 
logical models of hadrons. Of course, the leading order relations (i.e. those valid 
in the M — > oo limit) have to be corrected for symmetry breaking effects in order 
to connect them to properties of physical hadrons containing heavy quarks. These 
corrections arise from unequal quark masses and from terms suppressed by powers 
of 1/M, and are in the range of 10 - 20% for many instances involving b and c quarks. 
^ ■ Thus extracting the leading behaviour using heavy quark symmetries, and then es- 
timating the corrections using some phenomenological model, is a practical solution 
to cut down our ignorance in dealing with QCD and to make useful predictions for 
fitting experimental results concerning hadrons containing heavy quarks 2 . 

Of particular phenomenological interest are the weak decay form factors of 
hadrons containing a single heavy quark. Together with the experimentally ob- 
serverd weak decay matrix elements, these form factors determine various elements 
of the Cabbibo-Kobayashi-Maskawa quark mixing matrix. For example, a precise 
determination of the element V C b is possible provided we know the form factors of 
vector and axial currents between B and D mesons and B — > D semileptonic decay 
rates. The approach outlined above predicts that in the leading order all such form 
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factors can be reduced to two unknown functions, one for mesons (£) and another 
for baryons ((). These two we refer to as the Isgur-Wise functions. 

The symmetry properties of hadrons containing a single heavy quark are best 
expressed in the language of an effective field theory. In this formalism, factors of 
heavy quark mass are explicitly taken out by appropriately scaling all variables and 
the hadronic states are characterised by their four-velocities. In the M — > oo limit, 
the momentum carried by the QCD degrees of freedom (gluons and light quarks) 
that interact with the heavy quark is too small to alter v^, and so becomes a 
conserved quantum number. Of course, the large momentum transfer involved in a 
weak decay can change by a finite amount. £ and ( are thus functions of the only 
Lorentz invariant combination available, v ■ v', where v and v' are the heavy hadron 
velocities in the initial and final states respectively. (Note that v-v' is nothing but 
the Lorentz contraction factor 7 of special relativity, while v 2 = 1 = v' 2 .) 

In the rest frame of the heavy hadrons, the heavy quark just sits at one point, 
say the origin, while the light QCD degrees of freedom buzz around it. The situation 
corresponds to embedding a static colour impurity in the QCD vacuum. This picture 
provides explicit information about the functions £ and C in two specific geometries, 
v = v' and v = -v'. When v = v', the heavy quark may decay at some instance of time, 
but the surrounding light QCD degrees of freedom do not feel any change at all. 
The flavour independence of QCD thus allows a convenient choice for the absolute 
normalisation of the Isgur-Wise functions, £(v-v' = 1) = 1 = ((v-v' = 1). When v = -v', 
the actual physical process is heavy quark pair creation or annihilation. So we 
expect the Isgur-Wise functions to be singular at this kinematic point. Analytically 
continuing the Isgur-Wise functions to the full complex v-v' plane, a minimal scenario 
for their singularity structure is to have a pole at v-v' = -1, followed by a branch cut 
for v-v' < -1, and no other singularities in the rest of the complex plane. The physical 
situations of course correspond to \v ■ v'\ > 1 along the real axis. It is important to 
note that both these constraints on the Isgur-Wise functions remain unaffected by 
QCD renormalisation effects. An easy way to see this is to go to the temporal gauge 
A Q = 0, in which the static heavy quark loses all its QCD interactions and behaves 
like a free fermion. 

To go beyond these constraints and explicitly determine the Isgur-Wise functions 
requires knowledge of the dynamical behaviour of QCD. Formally the set up has 
several similarities to the Kondo problem of condensed matter physics, i.e. a static 
magnetic impurity embedded in a free electron conduction band. To name a few: 

• Both problems have a static impurity interacting with its polarisable surround- 
ings. Kondo problem has spin interaction, while heavy quark QCD has colour 
interaction. Only s-wave configurations feel the interaction in the leading order, 
and they are the cases of phenomenological importance. 

• The interaction is weak at short distances, but becomes strong at long distances 
completely screening the impurity. The change from asymptotic freedom to 
confinement is smooth without any intervening phase transition, and one can 
talk about renormalisation group (RG) flows leading from the weak coupling 
fixed point to the strong coupling one. 



• The unstable weak coupling fixed point is at zero coupling, with an essential sin- 
gularity and a logarithmically running coupling governing the scaling behaviour. 
The stable strong coupling fixed point is at infinite coupling, and is of a triv- 
ial nature. There is no inherent mass scale in either problem (the light quark 
masses can be set to zero in case of QCD). 

• The scale characterising the cross-over from weak to strong coupling is dynam- 
ically generated. This dimensional transmutation can be described as the ir- 
relevant operators determining the absolute scale of the theory. The Kondo 
temperature and A QC d are defined in terms of the couplings of an ultraviolet 
regulated theory. 

There are substantial differences between the two problems too: 
o In the Kondo problem, the electrons interacting with the magnetic impurity are 
free, with no interactions amongst themselves in the absence of the impurity. 
The electron energy spectrum extends continuously all the way down to E = 0. 
On the other hand, QCD is a non-trivial theory without any known solution 
even in the absence of a static impurity. Its states are discrete with a non-zero 
mass gap (ignoring the Goldstone bosons), 
o In the Kondo problem, the objects of study are the modifications of electronic 
properties caused by the static impurity. These can be studied in s-wave config- 
urations, while reducing the 3 + 1 dimensional problem to a 1 + 1 dimensional one. 
For heavy quarks in QCD, the objects of prime importance are not the changes 
in QCD excitation spectrum due to the heavy quarks, but the behaviour of the 
static colour sources themselves. This behaviour is too simple in the v = ±v' ge- 
ometries, as discussed above. Dynamical features of interest involve v ^ ±v', and 
these geometries cannot be reduced from 3 + 1 dimensions to 1 + 1 dimensions. § 
The modern solution to the Kondo problem is described in the language of 
conformal field theory. The differences listed above, however, make it unlikely that 
the problem of heavy quarks in QCD can be solved by the same techniques. On 
the contrary, the first solution to the Kondo problem was provided by Wilson in 
the framework of the renormalisation group 3 4 , and that is an approach which can 
definitely be applied to QCD. We thus turn our attention to the logic used by Wilson 
in solving the Kondo problem. 

The weak coupling and strong coupling expansions for the quantities of interest 
were known in the case of the Kondo problem. Wilson used RG to interpolate 
between the two, fixing the ratio of the dimensionful scale parameters appearing 
in the two solutions. This needed numerical RG integration keeping track of the 
low lying states of the problem. The desired electronic properties could then be 

§ It is possible to argue that arbitrary functions of v ■ v' can be Taylor expanded 
around the point v ■ v' = 1, and the derivatives can be evaluated in a purely s-wave 
geometry. We then have to calculate matrix elements of the original operators plus 
their higher dimensional descendents defined by insertions of covariant derivatives in 
the original operators (such towers of operators are familiar objects in the operator 
product expansion and conformal field theories). 



evaluated numerically. The problem of heavy quark QCD is simpler in a sense, be- 
cause it has operators which are RG invariant. Among these operators are partially 
conserverd vector and axial currents and their divergences — precisely the objects 
which are used to define the Isgur-Wise functions. The task is to implement the 
renormalisation group for QCD in such a way that the desired functions appear in 
RG invariant matrix elements, evaluate the matrix elements in the strong coupling 
limit, and then just read off the desired functions from the results by separating 
out the appropriate proportionality factors. 

This procedure can be followed in the lattice formulation of QCD 5 . The for- 
mulation is non-perturbative, yet at any non-zero value of the lattice spacing one 
can define a partially conserved vector current 6 . The RG invariant divergence of 
the vector current is just the difference of the quark mass operators. The partially 
conserved vector current is a non-local object on the lattice and its precise form 
depends on the details of the lattice action (it is a l-link long operator in the sim- 
plest lattice discretisations), but its divergence is a local operator whose structure 
can be chosen independent of the details of the lattice action. 

For concreteness, let us consider the vector current form factor for B — > D 
semileptonic decay and its limiting behaviour as the b and c quark masses go to 
infinity. 

< £>(j/)|c7„&|S(p) > = (p + p')Mq 2 ) + qj-(q 2 ) 

AI =T ^M B M D (v + v% £(v ■ v') . (1) 

Here q = p-p' is the momentum transfer, and v ■ v' = {M% + M% - q 2 )/2M B M D . The 
divergence of the above expression gives 

<D(p')\clb\B(p)> = (M B + M D )f + (q 2 ) + - q \ f /_(g 2 ) 



M B -M D n ' v/ v M B -M D ' 

M -^° y/M B M D (1 + v v 1 ) £(v ■ v') . (2) 

Note that the operators appearing in the matrix elements on the left hand side of 
these equations are composed of the QCD fields, and not the rescaled fields of the 
heavy quark effective theory. This feature is essential to keep the RG evolution of 
the matrix elements simple, e.g. both sides of the equations are RG invariant and £ 
extracted from the above form factors does not have anomalous scale dependence. 

These results are still in the continuum language. Following Bjorken, an upper 
bound on £'(l) is obtained by equating the inclusive sum of probabilities for decays 
into hadronic states to the probability for free quark transition. The matrix elements 
appearing in Eq. (1) give 7 ' 8 £'(1) < -1/4, while those appeaing in Eq. (2) yield a 
stronger constraint £'(1) < -1/2. The improvement is due to the fact that the right 
hand side has a kinematic factor of (1 + v.v') in the divergence equation instead of 
{v + i n the vector current case. 

Analogous expressions can be written down for form factors of baryons con- 
taining a single heavy quark. For instance, the vector current form factor in the 



semileptonic A fc — > A c decay is 



< A c (p')|c7 M 6|A fc (p) > = u(p')[/i7m - if2<r^<f + fsq^Hp) 



M- 



[u(v')^u(v)} C(w • «') , (3) 



where u(p') and u(p) represent the spinor wavefunctions of the spin-half baryons. 
The divergence of this expression yields 

n ' h " '"' < A c (p')\clb\A b (p) > = u( P ')[h + - q \ f h]u(p) 



M Ab ~M Ac " 1 M Ab -M A 

M -^° [u(v')u(v)] C(v ■ v') . (4) 

Unlike the meson case, there are no kinematic factors here accompanying the func- 
tion £ on the right hand side. Therefore, the sum rule analysis gives only a weak 
constraint £'(1) < 0. 

It is worthwhile to note that a degenerate heavy quark mass limit (i.e. m b = m c ) 
can be smoothly taken for both Eqs. (2) and (4) after cancelling out the mass 
factors on the left hand side. This limit makes the heavy quark flavour symmetry 
exact, and prevents anomalous dimension factors (e.g. functions that behave as 
powers of a s (m c )/a s (m h ) in leading order weak coupling perturbation theory) from 
appearing in the RG analysis of the form factors, thus simplifying the extraction of 
the Isgur-Wise functions. 

Now we put QCD on a Euclidean lattice using staggered fermions 9 . The residual 
chiral symmetries of this implementation protect the quark mass operator from un- 
wanted renormalisations. The formalism for performing strong coupling expansions 
in lattice QCD is well-known. To keep the matters simple, the expansion in l/g 2 is 
used often in conjunction with simultaneous expansions in 1/N C and l/d. We follow 
this approach, i.e. first obtain results in the limit of infinite number of colours and 
infinite number of space-time dimensions and then look at the corrections due to 
finite iV c and d. 

For our purpose, it suffices to look at the extreme strong coupling limit of the 
theory. It is described by a trivial fixed point of the RG evolution. This fixed point 
is reached by carrying to extreme the procedure of integrating out all the higher 
energy virtual states of the theory while lowering the cutoff. Only the lowest state, 
described by a perfectly screened delta-function wavefunction in position space, sur- 
vives in each quantum number sector. All the radial excitations orthogonal to the 
lowest state, corresponding to extended wavefunctions, disappear. The interactions 
amongst the surviving states are of course altered from their weak coupling be- 
haviour to compensate for the disappearance of the excited states. We evaluate the 
correlation functions in position space using the summation of hopping parameter 
expansion method 1041 to keep this intuitive picture clear. 

The form factors of interest are rather trivial to calculate in the strong coupling 
limit. The Feynman diagrams are more conveniently drawn in terms of colour 
singlet hadron lines rather than the original quark lines. Lattice artifacts show 
up in the formulae (e.g. sin(p) and cos(p) functions appear in momentum space 



propagators instead of p and p 2 ), but they are easy to keep track of. For N c — > oo, 
the 3-point correlation function between the external heavy hadrons and the quark 
mass operator corresponds to a tree level Y-shaped graph. For B — ► D matrix 
elements, the three mesons meeting at the vertex are the incoming B, the outgoing 
D and the scalar cb. The full correlation function is merely the product of the three 
meson propagators, with the sum in position space over all possible locations of the 
vertex producing the constraint of momentum conservation. 

In the d — > oo limit, the 3-point correlation function after putting the external 
hadrons on mass-shell is: 

G(%';t^ M ),%;0),%^-oo)) = ^ - B - - — {= — , (5) 

2 smh Ed 2smh_E B 1 + 2n b K c 2^ cos(g M ) 

where C's are the state normalisation constants, £'s are the energies and «'s are 
hopping parameters representing the renormalised quark masses (in the notation 
of Ref. 11). Upon amputating the external legs of the correlation function we 
get the desired matrix element, which is the last factor on the right hand side of 
Eq. (5). It can be converted to the continuum notation by identifying various lattice 
expressions. With 

M% -Ml = -(---) , m b -m c = U---) , M 2 C = — , (6) 

K \Kb K c J 2 \ Kb K c ) KbK c 

and the scalar meson mass M sc « M B + M D in the M — *■ oo limit, we obtain 

Mb-Md M b + Md 



< D(p')\clb\B(p) > = 



m b -m c 1 - q 2 /M% c 
M B - M D {M b + M D f ( , 1 



1 + 0(TF) • ( 7 ) 



m b -m c 2M B M D (l + v ■ v') \ M 

Comparing Eq. (7) and Eq. (2), and taking the degenerate heavy quark mass limit 
(i.e. M B = M D )i we identify 



Applying the same method to the quark mass form factor for the spin-half baryons, 
we obtain 

The difference between ^ and C arises entirely due to kinematic factors inherent in 
their definitions. 

This analysis shows that in the strong coupling limit the functional form of the 
Isgur-Wise functions is completely determined by the l/(M s 2 c -q 2 ) dependence, or 
scalar saturation, of the quark mass form factor. The remaining hadron mass factors 
just ensure that the Isgur-Wise functions are properly normalised, £(1) = 1 = C(l)- 
Indeed the fact that the normalisations turn out to be correct is a confirmation of 



the RG invariance of the form factor evaluation. The singularity at v ■ v' = -1 is as 
per expectations, but no branch cut shows up in this leading order result. 

It is straightforward to get rid of the d — > oo limit, and work directly in d = 4. The 
expressions for state normalisation constants, hadron masses, hopping parameters 
and correlation functions become different in terms of the bare parameters appearing 
in the lattice action 11 . The phenomenon of scalar saturation, however, is not altered 
at all and the results of Eqs. (8) and (9) remain valid. 

Scalar saturation no longer holds once processes suppressed by powers of 1/N C are 
taken into account. Such processes are represented by diagrams containing hadron 
loops. For example, to include the most dominant correction from dynamical light 
quark loops we have to evaluate the Feynman diagram with a heavy-light meson 
loop attached to three hadron legs, two corresponding to the external states and 
one corresponding to the quark mass operator. This process can be looked upon 
as the scalar cb decaying into two heavy-light mesons which in turn interact with 
the two external hadrons. Evaluation of such diagrams is quite involved, but we 
can estimate their magnitude by simple dimensional analysis. The virtual processes 
contained in the corrections suppressed by powers of 1/N C are strong interaction 
hadron vertices. The characteristic scale for these vertices (or decay widths) is 
Aqcd, which remains finite as M — > oo. As a result, the 1/N C corrections to scalar 
saturation are expected to be 0(T sc /M sc ) = 0(Aqcd/M). Such corrections do not 
contribute to the Isgur-Wise functions which are defined in the M — > oo limit. We 
thus argue that Eqs. (8) and (9) are correct even when N c is finite. A particular 
instance of 1/N C suppression is the branch cut for v ■ v' < -1, which we expect to be 
softer by a factor of A QC d/M compared to the pole at v ■ v' = -1. 

Now we can complete the task of inferring what form factors at the weak coupling 
fixed point of QCD could have evolved to the strong coupling expresssions obtained 
above. £ and C are dimensionless functions with fixed absolute normalisations at the 
no- recoil point, v and v' are not affected by RG evolution, and the operator involved 
in the calculation above (i.e. the local divergence of the partially conserved vector 
current in the limit of degenerate heavy quark masses) was carefully chosen to be 
RG invariant to avoid anomalous dimension corrections. We conclude that the 
results of Eqs. (8) and (9) are exact even for the weak coupling fixed point of QCD. 

We close with several comments regarding these results: 

(a) Since the quark mass operator and its RG evolution plays a crucial role in 
our analysis, staggered fermions were necessary in the lattice implementation. 
Wilson fermions 12 would not have been of much use, since in that case the quark 
mass operator undergoes additive renormalisation. 

(b) The strong coupling limit does not possess all the properties of the weak cou- 
pling fixed point of QCD, since different quantities follow different RG evolution 
patterns. The appropriate choice of an RG invariant quantity is a crucial in- 
gredient in connecting the two limits. Here the choice of the vector current was 
essential; the axial current does not have as nice renormalisation properties on 
the lattice. The Isgur-Wise relations amongst the form factors of the vector and 
the axial currents, which hold in the weak coupling limit, are not expected to 



hold in the strong coupling limit. This is not a disaster. As long as one of the 
form factors can be determined by connecting it from the strong coupling to 
the weak coupling limit, the rest can be fixed in the weak coupling limit by the 
usual continuum manipulations. 

Scalar saturation of the quark mass form factor is a simple consequence of strong 
coupling and iV c — > oo limits. Explicit lattice formulation is not necessary to infer 
this behaviour. One can anticipate it just on the basis of perfect screening 
between quark and antiquark, removal of all excited states in the process of RG 
evolution and suppression of couplings to multi-hadron virtual states. 

The series expansions in 1/g 2 and 1/M have non-zero radii of convergence, al- 
though the one in l/N c does not. This is sufficient to avoid any problems while 
interchanging the order of limits. Anyhow, the procedure followed here is to first 
calculate the form factor in the strong coupling limit, then take the degenerate 
quark mass limit, then let M become large to extract the Isgur-Wise functions, 
and finally worry about 1/N C corrections. 

The light quark mass is kept finite and constant throughout. It is not at all 
necessary to take the chiral limit for the light quarks. The fact that the light 
quark mass remains a constant in the M — > oo limit is enough, for instance to 
justify the replacement M sc w M B + M D . 

In principle, the complete form factors, i.e. without the use of M — > oo limit, can 
be evaluated in the strong coupling limit. The subleading corrections suppressed 
by powers of 1/M, however, are not universal and not easy to connect from the 
strong coupling to the weak coupling limits. For instance, the terms suppressed 
by 1/N C have to be kept, and the results have to be converted to continuum 
language using light hadron masses, hadronic excitation energies and widths, 
and so on. It may be possible to keep track of all this in an elaborate numerical 
RG evolution scheme involving many low lying states, such as the one employed 
by Wilson to solve the Kondo problem. 

The leading Isgur-Wise functions (which have £'(1) = -1 and C'(l) = -I/ 2 ) do a 
reasonable job in fitting the experimental data for semi-leptonic B-decays 1314 . 
A better check needs estimates of symmetry breaking corrections (from unequal 
heavy quark masses and from terms suppressed by powers of 1/M) that have to be 
added to the functions £ and C extracted above. With the leading term already 
taken care of, these corrections can be found using phenomenological models 
without introducing too much uncertainty in the final results. Quantitative 
estimates of the corrections together with more precise experimental data would 



provide an accurate test of the results presented here. 
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